Kondo resonance line-shape of magnetic adatoms on decoupling layers 
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The zero-bias resonance in the dl /dV tunneling spectrum recorded using a scanning tunneling microscope 
above a spin- 1/2 magnetic adatom (such as Ti) adsorbed on a decoupling layer on metal surface can be accurately 
fitted using the universal spectral function of the Kondo impurity model both at zero field and at finite external 
magnetic field. Excellent agreement is found both for the asymptotic low-energy part and for the high-energy 
logarithmic tails of the Kondo resonance. For finite magnetic field, the nonlinear fitting procedure consists 
in repeatedly solving the impurity model for different Zeeman energies in order to obtain accurate spectral 
functions which are compared with the experimental dl /dV curves. The experimental results at zero field are 
sufficiently restraining to enable an unprecedented reliability in the determination of the Kondo temperature, 
while at finite fields the results are more ambiguous and two different interpretations are proposed. 

PACS numbers: 72.10.Fk, 72.15.Qm 



When an atom with non-zero spin and/or orbital mag- 
netic moment is adsorbed on an ultrathin insulating ("decou- 
pling") layer grown on a metal substrate, it will behave as a 
nearly ideal magnetic impurity with well-defined local mo- 
ment which is described by the Kondo model Hj-Q]. The 
Kondo model explains the effects of the weak coupling of the 
impurity moment with the itinerant electrons in the substrate. 
If the coupling is antiferromagnetic, the impurity moment is 
screened in a complex many-particle process (the Kondo ef- 
fect J3-01), the signature of which is a resonance in the im- 
purity spectral function which peaks in close vicinity of the 
Fermi level li UTlll . The theoretically predicted resonance 
line-shape is parabolic near its maximum 0, Q, has long 
tails which can be empirically well fitted by an inverse square 
root (Doniach-Sunjic) function 11314 131 . but which have been 
shown to really be logarithmic jl6lll7ll : the Kondo resonance 
cannot be well approximated by a simple Lorentzian reso- 
nance BQ. 

Recent tunneling spectroscopy experiments performed at 
very low temperatures provide Al/dV spectra of magnetic 
adatoms with excellent energy resolution and very low noise 
1 2 11 - 42311 . Numerical techniques for accurately computing 
impurity spectral functions at finite frequencies have also 
reached maturity | L3, 15, 24-3(3]. We are thus presently able 
to directly compare experimental spectra with the theory. This 
is interesting for several reason. First, it is a non-trivial test 
of the adequacy of simplified quantum impurity models for 
describing the low-energy behavior of single magnetic adsor- 
bates. Second, it provides an unbiased approach for extracting 
the parameters of the impurity model. Third, it confirms our 
understanding of the logarithmic dynamical behavior of quan- 
tum impurities in the intermediate frequency regime (above 
the Kondo scale, but below the atomic scale). Our goal in 
this work is to fit the experimental tunneling spectra for a Ti 
adatom on the Cu 2 N/Cu(100) surface, Fig. 4 in Ref.|23, 



us- 



ing Kondo resonance curves calculated using the numerical 



which includes not only the Kondo resonance peak but also 
its long-tailed flanks. 

The experimental Kondo-resonance spectrum for Ti shows 
a sizeable asymmetry between positive and negative bias volt- 
ages and the resonance peak is slightly displaced from zero 
bias. These two features indicate possible presence of addi- 
tional potential scattering in the problem (in addition to the 
exchange scattering, which is the crucial element of the Kondo 
model). It is, however, also possible to account for the asym- 
metry using a model which takes into account quantum inter- 
ference between a narrow and a broad scattering channel (sim- 
ilar to the Fano formula used for fitting the spectra of magnetic 
impurities adsorbed directly on a metal surface Jiol [l8 . 321. 
but based on the correct impurity spectral function rather than 
on an approximation by a Lorentzian peak | i_8,|2(i), while the 
displacement of the peak might be explained by some small 
bias offset Vq for experimental reasons (such as thermovolt- 
ages generated by temperature differences in the STM head 
or in the cryostat wiring). As we show in the following, it 
is indeed possible to obtain an excellent fit to the experimen- 
tal data with the universal Kondo-resonance line-shape of the 
symmetric Kondo impurity model; it is not necessary to in- 
clude any potential scattering in the impurity model itself. 

The fitting is performed with the phenomenological ansatz 
function 

(dI/dV)(V) = a + b [(1 - q 2 )lmG{eV) + 2qReG(eV)] , 

(1) 

where the impurity Green's function is chosen to be 



G{ui) = G K ondo[(w - w )/A 



HWHMj 



(2) 



renormalization group (NRG) technique 11201. 13111 . The fitting 
window is the bias voltage window from -lOmV to lOmV, 



The ansatz function (Q~|) is chosen so as to reproduce the Fano 
formula for the case of a pure potential scatterer. The bias 
voltage is expressed in energy units as u> = eV, where e 
is the elementary charge; ljq = eVb with the voltage off- 
set Vo. Parameter a accounts for the possible featureless 
background conduction, parameters b and Ajjwhm set the 
height and width of the resonance, respectively, and q is a 



2 



Fano-model-like asymmetry parameter such that q = de- 
scribes a symmetric peak, while q ^ leads to an asymme- 
try. Finally, GKondo is the universal impurity Green's func- 
tion for the Kondo model, as computed using a highly accu- 
rate NRG calculation. It is normalized and rescaled so that 
in zero magnetic field A(lj = 0) = (— l/7r)Im GKondo = 
0) = 1 and A(u> = 1) = 1/2. We show in the following 
that the five-parameter ansatz in Eq. (Q3 is sufficient to ob- 
tain nearly perfect agreement between the experimental re- 
sults and the theoretical spectral function at zero magnetic 
field throughout an extended energy range including the res- 
onance peak and its tails. We perform nonlinear fitting using 
the Levenberg-Marquardt method to obtain the parameter set 
(a, b, Ahwhm, ?! wo). To quantify the "goodness of fit" we 
calculate squared residuals a 2 for all data points. 

We first discuss the Kondo resonance in the absence of the 
external magnetic field, B = 0. The comparison of the theo- 
retical function, Eq. (Q~|i, with the experimental data is shown 
in Fig. [T]and the extracted parameters are tabulated in TableQ] 
(first line). Some comments are in order. First, the asymme- 
try parameter q is non-zero, but small. It is not necessary to 
add some polynomial background to the model in order to de- 
scribe the asymmetry. A linear background hardly affects the 
goodness of fit a 2 . A quadratic background somewhat im- 
proves agreement in the tails, yet - importantly - the value of 
q as extracted using such an extended fit function changes by 
only 3%. Thus we conclude that the Fano interference is the 
correct physical interpretation of the observed asymmetry. To 
further support this claim, additional calculations have been 
performed for the Kondo model with potential scattering, the 
Kondo model with conduction-band density of states with a 
finite slope, as well as for the asymmetric single-impurity An- 
derson model; the results indicate that it is not possible to re- 
produce the observed asymmetry without including the Fano- 
like interference in the ansatz function, except for rather ex- 
treme (and thus unphysical) values of the parameters in those 
models. 

The second comment concerns the parameter Ahwhm- Ac- 
curate NRG calculations indicate that the relation between the 
HWHM of the spectral function for spin- 1/2 Kondo model at 
the particle-hole symmetric point and the Kondo temperature 
(as defined by Wilson) is Q30Q 



Ahwhm /TV, ~ 3.7. 
We thus find ksTj^.w = 0.176 meV or 



T K . 



W 



2.04 K. 



(3) 



(4) 



Away from the particle-hole symmetric point the ratio 
Ahwhm /Tk,w is different (it depends on the quasi-particle 
scattering phase shift) [f3oh - It is also worth mentioning that 
there are several different definitions of the Kondo tempera- 
ture for the 5 = 1/2 Kondo model which are in common 
use. The Wilson's definition Tk.w is often used in NRG, per- 
turbation theory, and Bethe Ansatz studies, another definition 
T^?' is used in local Fermi liquid works, and a third defini- 
tion Tk,h nas been used in the works of Hamann, Nagaoka 



and Suhl. They are related through Tk,w 
T K H = 2.2T K w . We thus find 



(0) 



4.94 K, 



K 

T KH = 4.49 K. 



0.4128T^ 0) and 



(5) 



Some care is needed when comparing works where the Kondo 
temperature is defined in non-equivalent ways. 

The final third comment concerns the thermal broadening. 
The calculation has been performed at T = 0, while the ex- 
periment 12311 is performed at T = 500 mK, which is smaller 
than Tk,w on ly by a factor of four. It is known, however, that 
the main effect of finite temperature in the T < Tk range is 
to reduce the height on the resonance peak and only slightly 
increase the peak width; the logarithmic tails are not affected 
significantly. This indicates that the finite-temperature effects 
are small and that the reliability of the extracted parameters is 
not reduced. The NRG calculation can be performed at finite 
temperature, but this does not improve the quality of the fit 
in the present case. Nevertheless, it would be interesting to 
repeat the experiment at even lower temperature in an attempt 
to achieve the true asymptotic zero-temperature limit and to 
reduce the noise even further. 



• experimental data 
- fit 
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Figure 1 : (Color online) Theoretical function (light line, red online) 
and experimental data (black symbols) for the Kondo resonance in 
the dl /AV spectral function for a Ti magnetic impurity atom ad- 
sorbed on the CU2N layer on Cu(100) surface. The theoretical func- 
tion is given by Eq. ifTJ, the parameters are explained in the main text. 
The experimental results have been extracted from Fig. 4 in Ref.l23l 

For reference, in Fig. [2] we plot a number of spectral line- 
shapes for Fano-like interference processes for cases where 
the narrow resonant scattering channel is described either by 
a Lorentzian curve or by a Kondo resonance curve (see also 
Refs. 2(3, 33). The width of the Lorentzian has been fixed 



so that in the asymptotic low-energy region both curves over- 
lap. As a general rule, at the same value of the asymmetry q 
the "Fano-Kondo" line-shapes have longer tails and smoother 
variation compared to their "Fano-Lorentzian" counterparts. 
In particular, it should be noted that the experimental dl jdV 
curve in Fig. [TJ cannot be well described using a Fano- 
Lorentzian line-shape. 
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Parameter a b q Ahwhm IjJo q^bB/Tk 

B = OT 0.0239 ± 0.0008 -0.00472 ± 0.00002 0.058 ± 0.003 (0.651 ± 0.008) meV (-0.116 ± 0.003) meV (*) 

B = 7T(A) 0.3112 ±0.001 -0.00417 ± 0.00005 0.021 ± 0.007 (0.651 ± 0.008) meV (*) (-0.16 ± 0.01) meV 5.14 

B = 7T(B) 0.3112 ±0.0005 -0.0052 ± 0.0002 0.029 ± 0.003 (0.409 ± 0.002) meV (-0.133 ± 0.003) meV 7.36 



Table I: Full set of fit parameters for the Kondo resonance peaks shown in Figs.[T]and[3] The error estimates indicated include only the standard 
deviation from the fitting procedure. The asterisk indicates a parameter whose value is fixed in the minimization. 




Figure 2: (Color online) Comparison of the resonance line-shapes 
arising from Fano-like quantum interference in processes where the 
resonant channel is described by a Lorentzian (dashed line, red on- 
line) or by the Kondo resonance curve (full line, black online). 



We now turn to the dl /dV spectra in the presence of the 
magnetic field, where the splitting of the Kondo resonance is 
observed I34l - l38ll . We need to establish the value of a new pa- 
rameter, the Zeeman energy b = g/igB. We note that while B 
is known from the experiment, the g-factor is not and needs to 
be dermined by the fitting procedure that we now perform. 
This is a very non-trivial task, since the Kondo resonance 
splitting is not linear as a function of the Zeeman energy lf34l — 
3a 39ll . For each value of the ratio b/Tx a different univer- 
sal (split) Kondo resonance curve needs to be computed l42ll . 
In minimizing the parameters in the model function, Eq. (Q]), 
we first fix Ahwhm: this choice is equivalent to assuming 
that the magnetic field does not change the effective exchange 
coupling Jk of the impurity spin, which is a reasonable as- 
sumption. The resulting fit is shown in Fig[3](top panel), the 
parameters are tabulated in Table U (second line, A), and the 
variation of the residual error as a function of b/Tx is plotted 



in Fig. HI The agreement is fairly good, but not outstanding: 
the amplitude of the two resonance peaks is clearly underes- 
timated. From the known field B = 7K and the extracted 
b/Tx ratio, we determine the value of the g-factor: 



2.2. 



(6) 



This value is in the same range as the g-factors for other mag- 
netic adatoms on the same surface 112 111 . 
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Figure 3: (Color online) Theoretical function (light function, red 
online) and experimental data (black symbols) for the Kondo reso- 
nance splitting in the external magnetic field B = 7 T. Upper panel: 
fixed parameter Ahwhm- Lower panel: unconstrained Ahwhm- 
The magnetic field strength is 7 B. The experimental data have been 
extracted from Fig. 4 in Ref.l23l. 

The fit may be improved by relaxing the constraint that the 
effective Kondo exchange coupling Jk does not vary with the 
magnetic field. The reason for doing so is the observation that 
the S — 1/2 behavior at low energy scales in Ti adsorbates 
originates from the freezing out of the S = 1 degree of free- 
dom on some higher energy scales; indeed, a Ti atom in free 
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Figure 4: (Color online) "Goodness of fit" a (sum of the residual 
squared errors) for the field-splitting of the Kondo resonance assum- 
ing constant or unconstrained parameter Ahwhm as a function of 
the Zeeman energy to Kondo temperature ratio. 



space has two 3d electrons and is in a S 1 = 1, L = 3 state 



|23|,|40[]. Upon adsorption on the surface, the local orbital 
moment is strongly quenched and the spin moment is also 
reduced. The presence of the external magnetic field might 
affect how the original impurity degrees of freedom are com- 
bined into an effective S — 1/2 object and how this effec- 
tive spin is coupled with the substrate electron 114 111 . Since 
the Kondo temperature depends exponentially on Jk, i.e., 
Tjc oc exp(— 1/pJk), where p is the density of states in the 
substrate, even a small variation of Jk might induce a sizable 
change in Tk- We thus performed a second fitting calculation 
where the parameter Ahwhm was allowed to change. The 
result of this fit is shown in Fig. [3] bottom panel, and the pa- 
rameters are listed in TableU second line (B). The agreement 
improves significantly, see also the residual errors plotted in 
Fig. [4] The Kondo temperature determined by this approach is 
smaller by 40% compared to Tk at zero field, and the g-factor 
is found to be g = 3.2. While the purported change in Tk in 
the field is quite large and the value of g is larger than typical 
values for other adatoms, it is not clear whether the improve- 
ment of the fit can be solely ascribed to having a further free 
parameter in the ansatz function (i.e., overfitting) or if there 
are physical grounds for reduced Jk in the effective model. 
This ambiguity calls for further systematic experimental stud- 
ies. 

It has been shown that the zero-bias peaks in S = 1/2 
impurities on decoupling layers can be accurately described 
by the universal Kondo spectral function, confirming the log- 
arithmic frequency-dependence in the resonance tails. The 
Kondo temperature has been extracted in a reliable way. Due 
to non-linear behavior of the magnetic -field splitting of the 
Kondo resonance, two interpretations of the peak structure 
have been proposed which differ in the values of the exchange 
coupling Jk and the g-factor. 

I acknowledge the support of the Slovenian Research 
Agency (ARRS) under Grant No. Zl-2058 and Program Pl- 
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